The Role of Monopoles for Color Confinement 



H. Ichie *, A. Tanaka and H. Suganuma a 

a Research Center for Nuclear Physics (RCNP), Osaka University 
Mihogaoka 10-1, Ibaraki, Osaka 567, Japan 

We study the role of the monopole for color confinement by using the monopole current system. For the self- 
energy of the monopole current less than ln(2d — 1), long and complicated monopole world-lines appear and the 
Wilson loop obeys the area law, and therefore the monopole current system almost reproduces essential features of 
confinement properties in the long-distance physics. In the short-distance physics, however, the monopole-current 
theory would become nonlocal due to the monopole size effect. This monopole size would provide a critical scale 
of QCD in terms of the dual Higgs mechanism. 



1. Kosterlitz-Thouless Type Transition in 
Monopole Current Dynamics 

Color confinement can be interpreted using 
the dual version of the superconductivity ]l[-Q. 
In this picture, condensation of color magnetic 
monopoles is the key concept and leads to squeez- 
ing of the color-electric flux between quarks 
through the dual Meissner effect. QCD is reduced 
into an abelian gauge theory with monopoles in 
the 't Hooft abelian gauge[Q]-||. As for the max- 
imally abelian (MA) gauge, recent lattice stud- 
ies suggest that only abelian gauge fields with 
monopoles are essential for description of nonper- 
turbative phenomena in the infrared region||- 

We study the role of monopoles for confine- 
ment using a monopole action. In the con- 
finement phase, only short-range interactions re- 
main among monopoles |9| due to the screening 
effect |l0| by the dual Higgs mechanism, and 
therefore the infrared monopole-current dynam- 
ics would be described by a simple action || 
[qJ, S , =cv lat J] SAl fc{f t (s) 2 in the lattice formal- 
ism. Here, fcjf t (s)=e/(47r) • a 3 k, 1 (s) &Z and 
a lat =a/2-(47r/ea) 2 denote the magnetic current 
and its self-energy on the lattice, respectively. 

The monopole current system is generated with 
the current conservation condition, d^k™* = 0. 
Here, we update the link of the monopole current 
using the Metropolis method. 

We show in Fig.l the monopole density 



PM=jy J2 S ii I^Ji at ( s ) wrtn the lattice volume 
V=8 4 and the clustering parameter n = 
I Li) 2 , where Li denotes the loop length 
of the i-th monopole cluster. With increasing 
a, the monopole density pu is reduced gradu- 
ally for a~a c . The clustering parameter is dras- 
tically changed at a c . Quantitatively, the critical 
value of monopole self-energy, a c ~1.85, is close 
to ln7~1.95, which is derived from the entropy 
factor (2d — \) L = 7 L in the 4-dimensional self- 
avoiding random walk |]ll| . Such a transition is 
quite similar to the Kosterlitz-Thouless type tran- 
sition in 1+2-dimensional superconductors. Vor- 
tex condensation at high temperature in the su- 
perconductor corresponds to monopole condensa- 
tion in the strong coupling region of QCD. 

2. Dual Field Formalism and Role of 
Monopoles for Confinement 

In this section, we study confinement proper- 
ties in the monopole current system with k^O 
and .7^=0, which is the dual version of QED. 
To estimate the abelian Wilson loop, we intro- 
duce the dual gauge field B^x), since the gauge 
field A^x) inevitably includes the singularity 
as the Dirac string in the presence of magnetic 
monopoles. Using the relations d^* F^ v — k v and 
*F llv =d IJi B v —d v B lx , the dual gauge field B^ can 
be described from the monopole current k^, 
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in the dual Landau gauge d^B^ — 0. 
The Wilson loop is obtained as (W) = 
(expiie^e^ap f da^(d a Bf3 - dpBa)}). In esti- 
mating the integral in Eq.(]l|) numerically, we 
should take a finer mesh than original lattice 
spacing a to extract B^ correctly. Otherwise, 
the unexpected electric currents appear as the 
lattice artifact. However, too finer mesh is not 
necessary, because the original current fcjf'fs) in- 
cludes the error in the order of a. Numeri- 
cal analyses show that the calculation with c < 
a/2 is good enough for the estimation of B^, 
and hence we take c = a/2 hereafter. On lat- 
tices, the dual gauge field 6>^ ual = aeB^/2 is de- 
fined in the dual Landau gauge; #^ ual ( s + /i) = 
2nJ2 s > A( s - S ')^ at 0') using the lattice Coulomb 
propagator A(s). The dual version of the abclian 
field strength 6»^ al = ea 2 *F^/2 is defined by 
6^ al (s) = <9; t 6^ ual (s) - diesis), where d' de- 
notes the backward derivative. 

As shown in Fig. 2, the Wilson loop (W) ex- 
hibits the area law, which leads to the linear con- 
finement potential. In Fig. 3, we show the string 
tension cr lat = a phys a 2 in the lattice unit. For 
a$a c , the string tension is finite due to the long 
monopolc currents, the current simulations for 
(X^a c indicate that tr lat ~0 in the system without 
monopole currents. On the other hand, lattice 
QCD simulations show that the long monopole 
currents cover the whole system and string ten- 
sion is finite in the confinement phase even for the 
large j3 as (3 = 2.3~2.5 on 16 4 lattices. There- 
fore, the monopole self-energy a extracted from 
QCD must be smaller than a c = \n(2d — 1) = ln7 
in the confinement phase. Here, to extract the 
monopole action, one should use the monopole 
part in QCD, where only monopole currents ex- 
ist. Since the electric current and the short-range 
Coulomb force cannot be expressed only by the 
monopole currents, it is dangerous to compare 
the monopole action with the full SU(2) QCD (or 
the abelian part); the former includes monopole 
currents only, while the latter includes also elec- 
tric currents and Coulomb force. In any case, 
the monopole action is to be extracted from the 
monopole part in QCD and the monopole self- 
energy a should be smaller than a c —lnl in the 



confinement phase. 

3. Monopole Size Effect and Critical Scale 
in QCD 

The aim of our study is to analyze the theoret- 
ical structure of QCD in terms of the monopole 
degrees of freedom. Our studies indicate that 
the multi-monopole system has a similar struc- 
ture to QCD in the MA gauge. However, the 
monopole in abelian-projected QCD has the un- 
derlying structure reflecting QCD. For instance, 
QCD-monopoles would be as extended objects like 
hadrons H , because it is not an elementary par- 
ticle in QCD, but a collective mode composed by 
gluons. We compare the monopole current system 
with QCD ||] especially in terms of the self-energy 
and the size of monopolesJl2[. 

In the static frame of the monopole with an 
intrinsic radius R, the spherical magnetic field is 
created around the monopole as H(r) = xj&r = 
, \ * for r > R and H(r) = , for r<R, with 
QCD running gauge coupling e(r)fll||. 

For a>R, the electromagnetic energy observed 
on the lattice around a monopole is roughly es- 
timated as Af(a)~f a °°d 3 xiH(r) 2 ~j5^, which 
largely depends on the lattice mesh a. Hence, 
the monopole contribution to the lattice ac- 
tion reads S=M(a)a-L, where L denotes the 
monopole current length measured in the lat- 
tice unit a. Therefore, M(a) is related to the 
monopole-current coupling a lat and f3=2N c /e 2 , 
a lat ~M(a)a~Jf^=Z-0 SXJ{2) . By the screening 

effect, this relation becomes a lat ~j-/3 SU (2) @, 
which is consistent with the numerical result as 
shown in Fig. 3. Thus, as long as the mesh is 
large as a>R, the lattice monopole action would 
not need modification by the monopole size effect, 
and a lat is proportional to (3. 

On the other hand, there exists a critical cou- 
pling a c ~ln7 in the current dynamics. Above a c , 
the lattice current action provides no monopole 
condensation and no confinement, while (3~^oo 
can be taken in the original QCD keeping the 
confinement property shown in Fig. 3. Such a 
discrepancy between (3 and a can be naturally 
interpreted by introducing the monopole size 
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effect, which shouid lead the nonlinear corre- 
spondence between a and (3. For a<R, the 
extended structure of the monopole can be ob- 
served on the lattice, and the monopole creates 
the electromagnetic energy M(a) ~ — e -i^ R 
on the lattice, which is almost a-independent. Ac- 
cordingly, the lattice action should be changed as 

S=M(i?)aE (M n (f - 1* - A) 

because of the self-avoidness originating from the 
monopole size R. Thus, the monopole-current 
theory becomes nonlocal in the UV region. 

The monopole size R provides a critical scale 
for the description of QCD in terms of the dual 
Higgs mechanism. In the infrared region as a>R, 
QCD can be approximated as a local monopole- 
current action ||], and the QCD vacuum can be 
regarded as the dual superconductor. On the 
other hand, in the ultraviolet region as a\R, the 
monopole theory becomes nonlocal and compli- 
cated due to the monopole size effect, and the 
perturbative QCD would be applicable instead. 
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for useful discussions. H.I. is supported by Re- 
search Fellowships of the Japan Society for the 
Promotion of Science for Young Scientists. 



Figure Captions 

Fig.l The monopole density pm and the cluster- 
ing parameter rj as the function of a. 
Fig. 2 The Wilson loop (W(I x J)} on the lattice 
with the mesh c = a/2 for a — 1.7, 1.8, 1.9. 
Fig.3 The string tension er lat = a phys a 2 as the 
function of a in the multi-monopole system. The 
dotted line denotes the Creutz ratio in the lattice 
QCD with (3 = ±a. 
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